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Abstract. Observations of the long-duration 
Gamma-Ray Burst GRB 080916C by the Fermi 
Gamma-ray Burst Monitor and Large Area Tele- 
scope show that it has a single spectral form from 
8 keV to 13.2 GeV. The E > 100 MeV emission 
was « 5 s later than the E <, 1 MeV emission 
and lasted much longer even after photons with 
E < 100 MeV became undetectable. The redshift 
from GROND of z w 4.35 means that this GRB 
has the largest reported apparent isotropic 7-ray 
energy release, E iso w 8.8 x 10 54 ergs. It also sets 
a stringent lower limit on the GRB outflow Lorentz 
factor, r ,nin ~ 890, and limits the quantum gravity 
mass scale, M QG > 1.3 x 10 18 GeV/c 2 . 
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burst, Radiation mechanisms, Gamma-ray telescopes 
and instrumentation 



I. Introduction 

Gamma-ray bursts (GRBs) are the most luminous 
explosions in the universe and are leading candidates for 
the origin of ultra high-energy cosmic rays (UHECRs). 
Prompt emission from GRBs from ^10 keV to ~1 
to 5 MeV has usually been detected but, occasionally 
photons above 100 MeV have been detected by the En- 
ergetic Gamma-Ray Experiment Telescope (EGRET)!]]], 
and more recently by Astro-rivelatore Gamma a Immag- 
ini LEggero (AGILE) [2|. Observations of gamma rays 
with energies > 100 MeV are particularly prescriptive 
because they constrain the source environment and help 
understand the underlying energy source. Although there 
have been observations of photons above 100 MeVQ, 
1 4 1 , 1 5 1 , it has not been possible to distinguish competing 
interpretations of the emission[6|, J7), 0. The Fermi 
Gamma-ray Space Telescope, launched on 1 1 June 2008, 
provides broad energy coverage and high GRB sensitiv- 
ities through the Gamma-ray Burst Monitor (GBM) and 
the Large Area Telescope (LAT)|9|. The GBM consists 
of 12 sodium iodide (Nal) detectors which cover the 
energy band between 8 keV and 1 MeV, and two bismuth 
germanate (BGO) scintillators which are for the energy 
band between 150 keV and 40 MeV. The LAT is a 
pair conversion telescope with the energy coverage from 
20 MeV to more than 300 GeV. In this paper, we report 
detailed measurements of gamma-ray emission from the 
GRB 0809 16C detected by the GBM and LAT. 



II. Observations 

At 00:12:45.613542 UT (T ) on September 16 2008 
the GBM flight software triggered on GRB 0809 16C. 
The GRB produced large signals in 9 of the 12 Nal 
detectors and in one of the two BGO detectors. Analysis 
of the data on the ground localized the burst to a Right 
Ascension (RA) = 08 h 07 m 12 s , Declination (Dec) = 
-61°18'00"n0|, with an uncertainty of 2.8° at 68% 
confidence level (C.L.). At the time of the trigger, the 
GRB was located ~48° from the LAT boresight and 
on-ground analysis revealed a bright source consistent 
with the GRB location. Using the events collected during 
the first 66 s after T , within 20° around the GBM 
burst position, the LAT provided a localization of RA 
= 07 h 59 m 31 s , Dec. = -56°35'24"0~D with a statistical 
uncertainty of 0.09° at 68% C.L. (0.13° at 90% C.L.) 
and a systematic uncertainty smaller than ~0.1°. 

Follow-up X-ray and optical observations revealed 
a fading source at RA = 07 ,l 59 m 23.24 s , Dec. = 
-56°38'16.8" (±1.9" at 90% C.L.)Q2| by Swift/X- 
Ray Telescope (XRT) and RA = 07 h 59 m 23.32 s , Dec. = 
-56°38'18.0" (±0.5")C3, EH) by Gamma-Ray Burst 
Optical/Near-Infrared Detector (GROND), respectively, 
consistent with the LAT localization within the estimated 
uncertainties. GROND determined the redshift of this 
source to be z — 4.35 ± 0.15 lfT5ll . The afterglow was 
also observed in the near-infrared band by the Nagoya- 
SAAO 1.4 m telescope (IRSF)lT6l. 

The lightcurve of GRB 0809 16C, as observed with 
Fermi GBM and LAT, is shown in Fig. Q] The total 
number of LAT counts after background subtraction in 
the first 100 s after the trigger was > 3000. For most 
of the low-energy events, however, extracting reliable 
directional and energy information was not possible. 
After we applied standard selection cuts|9| for transient 
sources with energies greater than 100 MeV and di- 
rections compatible with the burst location, 145 events 
remained (panel 4), and 14 events had energies > 1 GeV. 

Because of the energy-dependent temporal structure 
of the lightcurve, we divided the lightcurve into five 
time intervals (a,b,c,d,e) delineated by the vertical lines 
(Fig. [T]). The GRB lightcurve at low energy has two 
bright peaks, one between and 3.6 s after the trigger 
(interval 'a'), and one between 3.6 and 7.7 s (interval 
'b'). The two peaks are distinct in the BGO lightcurve, 
but less so in the Nal. In the LAT detector the first peak 
is not significant though the lightcurve shows evidence 
of activity in time interval (a), mostly in events below 
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Fig. 1: Lightcurves for GRB 080916C observed by the 
GBM and the LAT, from lowest to highest energies. The 
energy ranges for the top two panels are chosen to avoid 
overlap. The top three panels represent the background- 
subtracted lightcurves for the Nal, the BGO and the LAT. 
The top panel shows the sum of the counts, in the 8- 
260 keV energy band, of two Nal detectors (3 and 4). 
The second is the corresponding plot for BGO detector 
0, between 260 keV and 5 MeV. The inset panels give 
views of the first 15 s from the trigger time. In all cases, 
the bin width is 0.5 s; the per-second counting rate is 
reported on the right for convenience. 



100 MeV. Above 100 MeV, peak (b) is prominent in the 
LAT lightcurve. Interval (c) coincides with the tail of the 
main pulse, and the last two intervals reflect temporal 
structure in the Nal lightcurve and have been chosen 
to provide enough statistics in the LAT energy band 
for spectral analysis. The highest energy photon was 
observed during interval (d): E = 13.22i° 5° GeV. Most 
of the emission in peak (b) shifts toward later times as 
the energy increases (inset). 



III. Spectral Analysis 

We performed simultaneous spectral fits of the GBM 
and LAT data for each of the five time bins described 
above and shown in Fig.Q] GBM Nal data from detectors 
3 and 4 were selected from 8 keV to 1.0 MeV, as 
well as BGO detector data from 0.26 to 40 MeV. 
LAT photons were selected using the "transient" event 
classiO for the energies from 100 MeV to 200 GeV 
This event class provides the largest effective area and 
highest background rates among the LAT standard event 
classes, which is appropriate for transient sources with 
small background integration time. This combination of 
the GBM and LAT data results in joint spectral fits using 
forward-folding techniques covering over six decades of 
energy. 

The spectra of all five time intervals are well fit by 
the empirical Band function ifTTl which smoothly joins 




Fig. 2: Fit parameters for the Band function, a, (3 and 
E pea k as a function of time. 



10° 

io- 
ht* 




NAI_03 + 
NAI_04 c 

BG0_00 
LAT « 




+++ 



+ + -+ -■ 



10 ! 



10' 

Energy (keV) 



10 s 



10 s 



Fig. 3: Count spectrum for Nal, BGO and LAT in 
time bin (d): the data points have ler error bars while 
upper-limits are 2a. The histograms show the number of 
counts obtained by folding the photon model through the 
instrument response models. The black histogram shows 
the count models obtained by folding the best-fit Band 
photon function through the detector responses. A fit 
which is a combination of the Band function and a high- 
energy power law is shown with the cyan histogram. The 
residuals are calculated for the simple Band function fit 
in the bottom panel. 



low- and high-energy power laws. Fig. [2] shows time 
evolutions of fit parameters, a, (3 and Speak- The first 
time interval, with a relative paucity of photons in the 
LAT, also has the most distinct spectral parameter values. 
The low-energy photon index a is larger (indicating 
harder emission) and the high-energy photon index (3 
is smaller (indicating softer emission) - consistent with 
the small number of LAT photons observed at this time. 
After the first interval there was no significant evolution 
in either a or f3. In contrast, -E pca k, the energy at which 
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the energy emission peaks in the sense of energy per 
photon energy decade, evolved from the first time bin 
to reach its highest value in the second time bin, then 
softened through the remainder of the GRB. The higher 
£-poak and overall intensity of interval (b), combined 
with the hard value of [3 that is characteristic of the 
later intervals, are the spectral characteristics that lead 
to the emission peaking in the LAT lightcurve (Fig. [T). 

We searched for deviations from the Band function, 
such as an additional component at high energiesQ. 
Three photons in the fourth time bin had energies above 
6 GeV. We tried modeling these high-energy photons 
with a power law as an additional high-energy spectral 
component. Compared to the null hypothesis that the 
data originated from a simple Band GRB function, 
adding the additional power-law component resulted in 
a probability of 1% that there was no additional spectral 
component for this time bin; with five time bins, this 
is not strong evidence for any additional component. 
Fig. [3] shows the spectrum of this interval with two fit 
models, a Band function and a combination of the Band 
function and an additional high-energy power law com- 
ponent. Our sensitivity to higher-energy photons may be 
reduced at z ~ 4.35 through absorption by Extragalactic 
Background Light (EBL). Because the effect of various 
EBL models ranges widely (T8), fl9l , we do not use 
EBL absorption effects in our estimation of significance, 
which gives conservative estimates. 

IV. Long Lived Emission 

Although the lightcurves shown in Fig. Q] indie ate that 
during interval (e) the spiky structures typical of prompt 
GRB emission appear to be dying out, the emission 
persisted in some of the GBM Nal detectors at a low 
level out to nearly To + 200 s. Since the excess above 
background in the 12 Nal detectors occurred in the ratios 
expected for the geometry of the detectors relative to 
the burst direction and this type of low-level, extended 
emission is a known phenomenon in at least some long 
GRBs[20|, we associate it with the GRB and fit the 
spectrum with a power-law index of —1.92 ± 0.21. 
Emission beyond To + 200 s fell below the threshold 
of the GBM detectors. 

Long lived emission in the LAT band was searched 
for by a unbinned maximum likelihood fit of a power- 
law spectral function for a point source at the GROND- 
determined burst location using the event class with min- 
imum backgrounds ("diffuse"). The fits are performed 
for two time intervals T + 100-200 s and T + 200- 
1400 s for events within 15° of the GROND localiza- 
tion coordinates. The upper bound was chosen because 
after To + 1400 s, the GRB off-axis angle increased 
from 50° to 62° resulting in decreased effective area. 
Contributions from instrumental, Galactic, extragalactic 
components were included in the fit, as well as the bright 
source Vela (which is located 13° from the GRB). Fits 
in both time intervals show the presence of significant 
flux. For the final time interval, T + 200 to T + 1400 s, 



the fit yields a flux of (6.4 ± 2.5) x 10~ 6 7 enr 2 s" 1 
for E > 100 MeV with a power-law photon index of 
— 2.8±0.5 at a significance of 5.6a. If the position of the 
point source is left free instead of fixed to the GROND 
localization, the fit yields a source position of RA = 
07 h 57 m 33 s , Dec. = -57°00'00" with an uncertainty of 
0.51° at 90% C.L. This location is 0.45° from, and in 
agreement with, the GROND GRB position. To solidify 
the association of this extended emission with the GRB, 
we performed the same source detection procedure for 
data from To — 900 s to To and no emission was 
observed. A search for emission beyond To + 1400 s 
was also fruitless. 

We therefore associate this long-lived component with 
the GRB and include it as a sixth and seventh time 
interval for comparison with the early-time emission 
(Fig. 3J. In the LAT data, a constantly declining high- 
energy flux with a power-law decay index of —1.2 ±0.2 
is seen throughout Tq + 1400 s (red points, Fig. 0J. 
On the other hand, the flux in the GBM band shows 
a slower decay initially and an apparent break in the 
lightcurve at ~ To + 55 s. The power-law decay indices 
are approximately —0.6 and —3.3 before and after the 
break, respectively. 
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Fig. 4: Fluxes (top panel) for the energy range 50 — 
300 keV (shown in blue open squares) and above 
100 MeV (red filled squares), and power-law index as a 
function of the time from To to To + 1400 s (bottom 
panel, LAT data only). Red points are obtained by 
spectral fits of the LAT-only data for all time intervals. 
Blue points are obtained with the Band functions listed 
in Table 1 for the first 5 intervals and a power-law fit 
with index -1.90 ± 0.05 for the 6th interval. 



V. Discussions 

The Fermi observations of GRB 0809 16C show that 
the event energy spectra up to ~ 100 s are consistent 
with a single model (Band function), suggesting that a 
single emission mechanism dominates. A non-thermal 
synchrotron emission is the favored emission mechanism 
at keV to MeV energies[21|, and can indeed reach 
- 30(r/1000) [5.35/(1 + z)} GeVE). However, it 
should be accompanied by a synchrotron self-Compton 
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(SSC) spectral component produced from electrons that 
Compton upscatter their synchrotron photons to 7-ray 
energies potentially in the LAT energy band. The ap- 
parent absence of an SSC component indicates that the 
magnetic energy density is much higher than the electron 
energy density or that the SSC vF v spectrum peaks at 
3> 10 GeV and thus cannot be detected. It should be 
noted that sensitivity to a high-energy additional spectral 
component is reduced because EBL can absorb high- 
energy photons via pair-production interactions although 
the effect of the EBL cannot be estimated reliably due 
to large model dependence. 

The delayed onset of the GRB 0809 16C LAT pulse, 
which coincides with the rise of the second peak in 
the GBM light curve (see Fig. [D indicates that the 
two peaks may originate in spatially distinct regions. 
A 77 pair-production opacity effect may be ruled out 
since we do not observe any spectral cutoff in the 
combined GBM/LAT data. It is intriguing that long lived 
gamma-ray emission in the LAT band exhibits different 
temporal behaviors from those in the GBM band. In 
particular, our measurement indicates a temporal break 
in the GBM band in contrary to the continuous decay in 
the LAT band. Before our observations, a high-energy 
(100 MeV - GeV) tail was observed most clearly from 
GRB 940217(13] in observations by EGRET, which was 
not conclusive. The LAT high-energy tail may indicate 
cascades induced by ultrarelativistic ions accelerated in 
GRBs[8], or angle-dependent scattering effects ll23l . 

A measurement of redshift z « 4.35 by GRONDlfT31 
combined with Fermi measurements provides several 
kinematic constraints on this GRB. Between 10 keV 
and 10 GeV in the observer's frame, we measure a 
fluence / = 2.4 x 10~ 4 ergs cm~ 2 which gives at 
z ps 4.35 an apparent isotropic energy release for a 
standard cold dark matter cosmology with cosmological 
constant Qa = 0.73, f2 m = 0.27, and a Hubble's 
constant of 71 km s _1 Mpc -1 of E iso w 8.8 x 10 54 ergs. 
This is ~ 4.9 times the Solar rest energy, and therefore 
strongly suggests on energetic grounds, for any stellar 
mass progenitor, that the GRB outflow powering this 
emission occupied only a small fraction ( ;$ 10~ 2 ) of 
the total solid angle, and was collimated into a narrow 
jet. 

Given the intensity of observed photons, large bulk 
Lorentz factor T is required to avoid the attenuation 
of high-energy photons in a compact emission region 
expected from rapid variability ll24l . Using the Band 
function as the target radiation field and setting to 
unity the optical depth r 77 to 7-ray pair production 
attenuation of the highest-energy observed photon, we 
obtain T m i n w 608 ± 15 and 887 ± 21 in time bins d 
and b, respectively, for z « 4.35. Our limits are much 
higher than previous firm estimates of F m i n ks 100 
from GRB 990123 11251 . The high-energy photons and a 
large redshift can also provides an limit on the violation 
of Lorentz invariance expected from some quantum 
gravity models [(26). In the linear approximation, the 



difference in the arrival times At is proportional to 
the ratio of photon energy difference to the quantum 
gravity mass, AE/Mqq, and depends on the distance 
the photons traveled. The arrival time of the 13.2-GeV 
photon relative to To, t — 16.54 s, is a conservative 
upper limit on its At relative to ~ MeV photons, and 
implies a robust lower limit on the quantum gravity 
mass, Mqg > 1-30 x 10 18 GeV/c 2 , which is an order of 
magnitude higher than the previous limit obtained in this 
fashion[ 27 1 . This lower limit is only one order of magni- 
tude smaller than the Planck mass, 1.22 x 10 19 GeV/c 2 . 

VI. Acknowledgement 

The Fermi LAT Collaboration acknowledges support 
from a number of agencies and institutes for both 
development and the operation of the LAT as well 
as scientific data analysis. These include NASA and 
DOE in the United States, CEA/Irfu and IN2P3/CNRS 
in France, ASI and INFN in Italy, MEXT, KEK, and 
JAXA in Japan, and the K. A. Wallenberg Foundation, 
the Swedish Research Council and the National Space 
Board in Sweden. Additional support from INAF in Italy 
for science analysis during the operations phase is also 
gratefully acknowledged. The Fermi GBM Collaboration 
acknowledges the support of NASA in the United States 
and DRL in Germany, and thanks Lisa Gibby, Al English 
and Fred Kroeger. 

References 

[1] B. L. Dingus, Astrophys. Space Scl, vol. 231, p. 187, 1995. 
[2] A. Giuliani et al, Astron. Astrophys., vol. 491, p. L25, 2008. 
[3] K. Hurley et al, Nature, vol. 372,' p. 652, 1994. 
[4] D. N. Wren, D. L. Bertsch, and S. Ritz, Astrophys. J., vol. 574, 
p. L47, 2002. 

[5] M. M. Gonzalez, B. L. Dingus, Y. Kaneko, R. D. Preece, C. D. 
Dermer, and M. S. Briggs, Nature, vol. 424, p. 749, 2003. 

[6] J. Granot and D. Guetta, Astrophys. J., vol. 598, p. Lll, 2003. 

[7] J. Katz, Astrophys. J., vol. 432, p. L27, 1994. 

[8] C. D. Dermer and A. Atoyan, New J. Phys., vol. 8, p. 122, 2006. 

[9] W. B. Atwood et al, "The large area telescope on the GLAST 
mission," submitted to Astrophys. J., 2008. 
[10] A. Goldstein and A. van der Horst, GCN Circular 8245, 2008. 
[11] H. Tajima et al, GCN Circular 8246, 2008. 
[12] M. Perri et al, GCN Circular 8261, 2008. 
[13] C. Clemens et al, GCN Circular 8257, 2008. 
[14] C. Clemens et al, GCN Circular 8272, 2008. 
[15] J. Greiner et al, Astron. Astrophys., vol. 498, p. L89, 2009. 
[16] T. Nagayama et al, GCN Circular 8274, 2008. 
[17] D. Band et al, Astrophys. J., vol. 413, p. 281, 1993. 
[18] S. Razzaque, C. D. Dermer, and J. D. Finke, astro-ph 0807.4294, 
2008. 

[19] F. W. Stecker, M. A. Malkan, and S. T. Scully, Astrophys.]., vol. 

648, p. 774, 2006. 
[20] V. Connaughton, Astrophys. J., vol. 567, p. 1028, 2002. 
[21] P. Meszaros, Annu. Rev. Astron. Astrophys., vol. 40, p. 137, 2002. 
[22] A. A. Pe'er and E. Waxman, Astrophys. J., vol. 613, p. 448, 

2004. 

[23] X. Y. Wang, Z. Li, and P. Meszaros, Astrophys. J., vol. 641, p. 
L89, 2006. 

[24] J. H. Krolik and E. A. Pier, Astrophys. J., vol. 373, p. 277, 1991. 
[25] Y. Lithwick and R. Sari, Astrophys. J., vol. 555, p. 540, 2001. 
[26] G. Amelino-Camelia et al, Nature, vol. 393, p. 763, 1998. 
[27] J. Albert et al, Phys.Lett.B, vol. 668, p. 253, 2008. 



